Scaling properties of critical bubble of homogeneous nucleation in stretched fluid of 
square-gradient density-functional model with triple-parabolic free energy 
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The square-gradient density-functional model with triple-parabolic free energy is used to study 
homogeneous bubble nucleation in a stretched liquid to check the scaling rule for the work of 
formation of the critical bubble as a function of scaled undersaturation A/j,/ A/ispin, the difference 
in chemical potential A/x between the bulk undersaturated and saturated liquid divided by A/ispin 
between the liquid spinodal and saturated liquid. In contrast to our study, a similar density- 
functional study for a Lennard- Jones liquid by Shen and Debenedetti [J. Chem. Phys. 114, 4149 
(2001)] found that not only the work of formation but other various quantities related to the critical 
bubble show the scaling rule, however, we found virtually no scaling relationships in our model near 
the coexistence. Although some quantities show almost perfect scaling relations near the spinodal, 
the work of formation divided by the value deduced from the classical nucleation theory shows no 
scaling in this model even though it correctly vanishes at the spinodal. Furthermore, the critical 
bubble does not show any anomaly near the spinodal as predicted many years ago. In particular, 
our model does not show diverging interfacial width at the spinodal, which is due to the fact that 
compressibility remains finite until the spinodal is reached in our parabolic models. 

PACS numbers: 47.55.db, 64.60.qe, 82.60.Nh 



I. INTRODUCTION 

Homogeneous bubble nucleation (cavitation) has at- 
tracted much attention for many years from a fundamen- 
tal point of view as well as for technological interests. 
It occurs when a fluid is undersaturated and is held at 
a pressure lower than its coexisting vapor pressure at a 
given temperature. Then the fluid is said to be stretched. 
Homogeneous bubble nucleation occurs in the stretched 
liquid. Homogeneous nucleation including homogeneous 
bubble nucleation in the liquid^ as well as liquid droplet 
condensation^ or the crystallization^ from melts or glass 
in general is the basic mechanism to initiate the first- 
order phase transition, and has been studied for many 
years mainly using a very simple analytical theory called 
classical nucleation theory (CNT)ii^i^. 

Recently, due to the development of the density func- 
tional theorji^ the deviations from the CNT has been 
predicted using various model systems These pre- 
dictions have successfully explained the deviation of the 
experimental data from the CNT predictions^^. There- 
fore, a reliable yet handy empirical formula which can 
predict nucleation rates is highly desired. To this end, 
various empirical formulae for bubble and droplet nucle- 
ation in a fluidiiii^'i^iii^'i^ as well as crystal nucleation 
in meltsi^iii have been proposed. 

The central quantities of the nucleation phenomena is 
the nucleation rate J, which is the number of critical 
nuclei formed per unit time per unit volume. Usually it 
is written in Arrhenius form 



J — A exp 



W* \ 



(1) 



where ^ is a kinetic pre-exponential factor which is be- 



lieved to be weakly dependent on temperature T, /cb is 
the Boltzmann's constant, and W* is the reversible work 
of formation of the critical bubble. Therefore, the tem- 
perature dependence of the nucleation rate is controlled 
by the work of formation W* , which most of the theoret- 
ical works study. 

The density functional theory is a powerful tool to 
study the density profile and the work of formation W* 
of critical nucleus. For example, Oxtoby and EvanSii 
has studied the work of formation W* and the density 
profile of critical droplets and bubbles of the Yukawa 
fluid using the so-called non-local density functional the- 
ory. Zeng and Oxtobj*^ studied the critical nucleus of 
the Lennard-Jones fluid using the non-local density func- 
tional theory. The present author^ used the local square- 
gradient density functional theory^'^'^ combined with the 
double-parabolic free energy, and demonstrated that the 
results of Oxtoby and coworkers^ can be reproducible 
by the square-gradient density functional theory rather 
than fully non-local theory. A similar model with triple- 
parabolic free energy is used later by BarrettiS,. Very 
recently, Li and Wilemsky-^ has compared the results 
obtained from the accurate non-local density functional 
theory with that from the approximate square-gradient 
theory and found that the two results agree qualitatively 
well. 

Among various nucleation processes, bubble nucleation 
has attracted relatively little attention^i^i^i^S. Recently, 
Shen and Debenedetti^^ studied homogeneous bubble nu- 
cleation in a stretched Lennard-Jones fluid using the non- 
local density functional theory of Zeng and Oxtobyi^. 
They found that if scaled by appropriate quantities, 
various quantities of the critical bubble including the 
work of formation W* degenerate into universal curves 
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when they are plotted against the scaled undersaturation 
A/i/A/igpin. Their results justify various empirical scal- 
ing approacbiiii^iii^ii^ii^ to the nucleation rate. Similar 
density functional calculations have been conducted by 
Kusaka not only for the Lennard- Jones fluid but for the 
square-well fluid. They found that the scaling is mode- 
latedly successfull for those two fluid systems. Recently, 
Punnathanam and Corti^^'^^ studied the cavity forma- 
tion rather than the bubble formation and found a similar 
scaling rule. 

All those studies mentioned above are based on the nu- 
merical results for several specific models of simple fluids. 
However, since the nucleation is so general and ubiqui- 
tous phenomena not restricted to simple fluid that it is 
highly desirable to study the scaling rule in more general 
models which can be applied not only to the simple liq- 
uid but to complex fluids, or even to amorphous or liquid 
metals. 

In this paper, we will use a simple square-gradient den- 
sity functional theory with a triple-parabolic free energy 
proposed by Granasy and Oxtoby^- to study the various 
properties of the critical bubble of homogeneous bubble 
nucleation. We choose this model as it captures the most 
basic properties of nucleation, yet many physical quan- 
tities can be handled analytically. In this study we par- 
ticularly focus on the scaling rule of various quantities'^. 
In Section II we present a short review of the square- 
gradient density-functional model with triple-parabolic 
free energy^! to summarize the necessary formula. We 
also correct a few typographical errors in the original 
article^!. In Section III, we will present the numerical 
results and discuss the implication of the results in light 
of the scaling rule. Finally Section IV is devoted to the 
concluding summation. 



which consists of three parabola centered at the vapor 
density 0o, and at the free energy barrier (pi (spinodal), 
and at the liquid density (p2, which we call "vapor" , "spin- 
odal" and "liquid" part of the free energy. We have cho- 
sen this triple-parabolic modeP^ rather than a simpler 
double-parabolic model^ as the former seems more real- 
istic near the spinodal point. 

Although Granasy and Oxtoby^i set (po = and 02 = 1 
to simplify the various formulae, we will leave two quan- 
tities finite as we want to consider the vapor phase with 
finite density 4>o- The curvature of parabola Aq and A2 are 
related to the compressibility of vapor and liquid phases, 
and A^ is the free energy difference between the liquid 
and the vapor. Although A^ represents in fact the pres- 
sure difference AP as Eq. ^ is the grand potential of 
open system', we call A/i chemical potential to make the 
comparison to the previous wor h-'^^i^'i^^i^'^i^^i'^i^'i'^ eas- 
ier since AP is proportional to A/i through AP = p\Afi 
with pi being the number density of the liquid^i^. We 
use the terminology "over-saturation" when Ap is pos- 
itive and "undersaturation" when Afj, is negative. The 
stretched liquid in this study corresponds to the under- 
saturated liquid. Therefore we will be mainly concerned 
with the situation when A/i < 0. 

From the continuity of the free energy Aaj(0) at the 
boundary (j)A and 0b, they are given by 



Ao0o + |Ai|0i 
Ao + IAil ' 

A202 + |Al|01 

A2 + IA1I 



(4) 



In Eq. ([2]) ^0 is the barrier height which separates the 
vapor at 00 from the liquid at (/12, and is given by 



II. SQUARE-GRADIENT 
DENSITY-FUNCTIONAL MODEL 

A. Triple-Parabolic Free Energy 

In the square-gradient density-functional model of the 
fiuid^-^'^, the free energy of the inhomogeneous fluid, such 
as the critical bubble in the stretched liquid is written as 



W 



Auj{(P) + c{W(pf) d 



(2) 



This form of the free energy is also known as the Cahn- 
Hilliard mode l'^i'^ or the phase-fleld mode l'^i'^ . In the 
triple-parabola model of Granasy and OxtobySi, the local 
part of the free energy Alo is given by 



Aw(0) = 



f (< 


" 0O)' 











A/i 



02 — 00 



A/i -I- ujo 



*> < (pA 

(3) 



< 



Wo 



A/i- 



(01 - 02) + A/i 



A/i, 



(5) 



where 



P = 
q = 



AolAi 



Ao + IAil' 

A2IA1I 
A2 + IA1I' 



(6) 



are constants determined from the compressibility A^. By 
equating two formulae for loq in Eq. ([5]), we obtain the 
location of the free energy barrier 



902) + Vm(0o - 02) -2Ap{p-q) 



(7) 



which depends on the over-saturation A/i > for the 
0)S'ersaturated vapor and the undersaturation Ap, < 
for the stretched liquid. The vapor spinodal is defined 
when the metastable vapor phase at (j)o becomes unstable. 
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This is realized when = 4>o, which leads to the over- 
saturation for the gas spinodal 



A^gas ^ -{(j)2- 00)" 



(8) 



Similarly the liquid spinodal occurs when the undersat- 
uration is given by 



P 



A^spin = (02 - M 



(9) 



free energy when Afi/ Afigpin = 0.3 near the coexistence. 
However, as the undersaturation increases (| A/uj becomes 
large) and it approaches the liquid spinodal A/Ltgpin, the 
density profile consists of only the two parts correspond 
to the spinodal and the liquid parts. We use the terminol- 
ogy " CNT regime" for the former regime where the clas- 
sical nucleation theory (CNT) is expected to be qualita- 
tively correct, and "spinodal regime" for the latter where 
the spinodal nucleation^ i^^'^°i'^^ is expected to occur. 



for the stretched liquid. Since we are interested in the 
critical bubble formation in the stretched fluid, we will 
consider the region from coexistence A^ = to the liq- 
uid spinodal A/igpin < 0. In contrast to the previous 
models^'SSii^'^ where the compressibility diverges con- 
tinuously as the spinodal is approached, the compress- 
ibility remains finite until the spinodal point is reached 
in our triple-parabolic model as the curvature Aq and A2 
is fixed. 
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FIG. 1: (a) The triple-parabolic free energy from CNT regime 
near the coexistence to the spinodal regime near the liquid 
spinodal for the case-I of section III (Table |l]| . (b) The cor- 
responding critical bubble at the CNT regime (A/j./ A/ispin = 
0.3) and at the spinodal regime (0.7). 

In Fig. [Tfa) we show typical shapes of the triple- 
parabolic free energy Auj together with the correspond- 
ing critical bubble in Fig. [Hb). The radii and 
are the matching radius that satisfies (va) = (pA and 
'/'(^b) = '/'B- Since the free energy consists of three 
parabolas corresponding to the vapor, spinodal and liq- 
uid parts, the density profile of the critical bubble con- 
sists of three parts that correspond to the three parts of 



B. Density Profile of Critical Bubble 

1. CNT regime 

Density profile of spherically symmetric critical bub- 
ble can be obtained from the Euler-Lagrange equation 
SW/Scj) = 0, which leads to the differential equation 



r 



0, 



0,1,2 



(10) 



for the three parabolas in Eq. ((2]), where F^ = ■\/|Ai|/2c 
and $i(r) = 0(r) — 0i, and -|- sign is used for i — \ and 
— is used for i = 0, 2 for ±. These differential equations 
should be solved with boundary conditions: 
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02 
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by 



together with the condition for $0 at the origin given 



dr 



= 



(12) 



for the critical bubble. In this CNT regime near the two- 
phase coexistence, both the matching radius ta and tb 
exist. When the under saturation A/i(< 0) is further 
increased, the matching radius ta approaches zero and 
disappears. Then only the matching radius tb can exist. 
We have named this regime near the spinodal point the 
spinodal regime. 

The solution of this Euler-Lagrange equation in 
Eq. pl?|) for the critical bubble are given by 

$o(0 = $oA?'yicsch(ForA)sinh(For)/r 

$i(r) = csc(Fi (rA - rs)) 

x(-$iBrBsin(Fi(r-rA)) (13) 
+^\ATa sin (Fi (r - rB)))/r 

'^2{r) ~ <i>2B?'B Gxp (— F2r + F2rs) /r 

For the three parts i = 0, 1, 2 respectively. 
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Finally, the matching radii rA and rs are determined 
from the simultaneous equation 



(14) 



where only $i is the functions of both rA and rs- 
This simultaneous equation can be solved numerically 
using standard algorithms such as the Newton-Raphson 
method. 



2. Spinodal regime 

In this case rA becomes zero. Therefore the i — Q part 
of the free energy density in Eq. ^ and its solution $o 
disappears. Then the boundary condition Eqs. (jlip and 
(fT^ for $0 should be replaced by the same condition for 
$1. The solution for the Euler-Lagrange equation for $2 
is the same as Eq. (|13l) . but the one for $1 now reads 



$i(r) ^ ^lijrscsc (Eirs) sin (Eir) /r, (15) 

for the critical bubble in the spinodal regime. 

In this case, the matching radius rs is simply deter- 
mined from the equation 



(16) 



+ IA1I, (17) 



d'^2 d<^i 

dr ^^^^ dr 

which is explicitly written as 

re (A2riC0t {T ir b) - \\i\T 2) = M 

that does not depend on the understaturation A/x. 
Therefore, the matching radius rs is constant in the 
spinodal regime. Again this equation should be solved 
numerically. 

In Fig. [TJb), we showed the typical density profiles in 
the CNT regime and in the spinodal regime. The critical 
bubble is larger in the CNT regime than in the spin- 
odal regime. However, the density difference between 
the inside and the outside of the bubble becomes smaller 
in the spinodal regime than in the CNT regime. Cor- 
respondingly, the interfacial thickness looks diffuse^^ as 
the spinodal is approached. 



derived by Cahn and Hilliard^^. 

Inserting the density profile Eq. ^ in the CNT 
regime, we can calculate the integral in Eq. (fT8|) ana- 
lytically, and we obtain 



W ^Wq + Wi+W2 



(19) 



where 



w 47r 3 , „ , , (l-rorACOth(rorA))$oArA 



3 ■ ■ V 92-90 

, 27r|Ai|0 



^ 

01 - 00 

A/i- — + A^ + uja 



— {-'^lArA + ^IBTB 



(20) 



W2 



- Fi ($iAri -I- $isr|) cot (Fi (r^ - re)) 

- Fi ($iA + ^ib) rATBcac (Fi (r^ - rs))) 

2TT\2(j)2^2BrB 



F2 
^ 2 



(l + rBF2) 



are the contributions of the three parts (vapor, spinodal, 
liquid) of the free energy density in Eq. Note that 
A/x < for the critical bubble in the undersaturated 
liquid. 



2. Spinodal regime 

Since, we only have the solution <I>2 in Eq. and $1 
given by Eq. (fT5|) for the critical bubble in the spinodal 
regime, we have 

W = Wi+W2 (21) 
where W2 is given by Eq. (^0)1 but Wi is given by 
47r . ' 



Wi 



A^- — + A^ + 0^0 



+ 



3 ' V 92-90 
27r|Ai|0i$iBri3 



(1-Firs cot (Firs)) (22) 



which can also be derived by setting = in Eq. ((20|) . 



3. work of formation of classical critical bubble 



C. Work of Formation of Critical Bubble 

1. CNT regime 

Once we know the density profile of the bubble and 
cavity, it is straightforward to calculate the work of for- 
mation W* of the critical bubble. To this end, we can 
use the formula 

H- = ^°°4„=(A.-i*^)* (18) 



Here, we will summarize the standard formula for the 
work of formation of critical bubble from the classical 
nucleation theory (CNT)^-^-'^. Using the capillary ap- 
proximation, the work of formation of critical bubble is 
given by the sum of volume term and the surface term: 

VF = — -r^lA^I +47rr27, (23) 
o 

where the critical bubble is assumed to be a sphere with 
the classical radius r with the sharp interface of the 
macroscopic surface tension 7. By maximizing this work 
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of formation with respect to the bubble radius r, we ob- 
tain the critical radius r^i of the critical bubble and its 
work of formation W* as 



'"cl = 



27 



(24) 



where 7 should be calculated at the liquid-vapor coexis- 
tence when Afj, = 0: 



2V-C 

7o + 



00 

7i - 



V Au!d(j) 

72 



(25) 



which consists of three contributions from the three parts 
of the free energy density. Due to the simple parabolic 
free energy for Aw, Eq. (pS)) can be calculated analyti- 
cally, and the results are!^ - 
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FIG. 2: Three shapes of triple-parabolic free energy corre- 
sponding to the three sets of free energy parameter listed in 
Table U studied in this article at the two-phase coexistence. 



70 = -V2Aoc((?!)A - (/-o)^ 



71 



72 



— arcsin 
2 



(26) 



-\/l\2C {(f>2 - 



Equation ([24ll reduces to the standard expression of the 
CNT~»2, if we recover the original definition A/i in Eq. ^ 
and replace A/j, by pi A/i or by the pressure difference AP. 



TABLE I: Three sets of free energy parameters used in this 
work to check the universality of the scaling properties of the 
critical bubble. 



model 


c 


00 


02 


Ao 


Ai 


Aa 


case-I 


1.0 


0.1 


1.0 


1.0 


-0.5 


0.3 


case-II 


1.0 


0.1 


1.0 


1.0 


-1.0 


0.7 


case-III 


1.0 


0.1 


1.0 


1.0 


-2.0 


0.9 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In order to check the universality of the scaling rela- 
tions for the various properties of critical bubble found 
by Shen and Debenedetti^i, we use this square-gradient 
density-functional model with triple-parabolic free en- 
ergy and study the various properties of the critical bub- 
ble. We use several typical free energy parameters to 
check the universality. Fig. [2] shows the three shapes of 
triple-parabolic free energy used in our work. Since the 
vapor is more compressible than the liquid in general, we 
choose Aq = 1 and A2 < Aq. The three sets of the free 
energy parameters used are summarized in TablelH Since 
our model is specified only by the compressibilities Ai and 
the undersaturation A/i, we will not consider the temper- 
ature effect^i^i^ explicitly for it will come in through the 
temperature dependence and the interrelations between 
A.. 

Figure [3] shows the matching radius ta and rs as 
the functions of the scaled undersaturation A/i/A/igpin 
for casc-I. The matching radius va becomes zero as the 
undersaturation enters the spinodal regime, while the 
matching radius rs survives and it becomes constant 
and independent of the undersaturation in the spinodal 
regime as predicted from Eq. ([T7)) . We also show the clas- 
sical radius r^i of the classical critical bubble predicted 
from Eq. ^ of the CNT. 



In Fig. [3] we also show the equimolar dividing radius 
r£i determined from 

pro roa 

/ ((/)(r) - (?!)org) 47rr^rfr = / {(j)2 - (j){r)) inr"^ dr 

Jo Jrn 

(27) 

where 

(/•org = ^ 0) (28) 
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FIG. 3: The matching radii rA and rs as well as the clas- 
sical radius r^i from the CNT for the case-I. Also shown is 
the equimolar dividing radius ro- In the spinodal regime ta 
disappears, radii rs and ro become constant. 
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is the density at the origin of the bubble. We use this 
definition Eq. (P?]) of the dividing surface used by Shen 
and Debenedetti^i though it is shghtly different from the 
usual definition of the dividing surface where (/<org should 
be replaced by (jjQ. Since the density profile (f)(r) is given 
by the analytical formulae in Eqs. (fT3|) and (fT5|l . not only 
the density at the origin (poig in Eq. ((28|) but the integrals 
on the both side of Eq. (P7|) are given by the analytic 
formula. Then, the non-linear equation Eq. (p7|) for the 
equimolar dividing radius ro can be solved numerically. 
By a careful inspection of Eq. (|27|l in the spinodal regime, 
we notice that the dividing radius rjj is constant as the 
matching radius rs in the spinodal regime. 

From Fig. [3] we notice that even though the classi- 
cal radius is continuously decreasing as the under- 
saturation |A/i| increases, the matching radius dis- 
appears while the matching radius and the dividing 
radius rjj becomes constant in the spinodal regime up 
to the liquid spinodal A/igpin- It is commonly believed 
that the radius of the critical nucleus diverges as the 
spinodal is approached^^'^^. However, this conclusion is 
derived from some mathematical models and could be 
model-dependent. In fact, this divergence is related to 
the third derivative of the free energy in the classical 
work of Cahn-Hilliard^Si. Also, the conclusion of diverg- 
ing critical nucleus drawn by Unger-Klein^S is based on 
the third expansion of the free energy that contains up 
to This third derivative does not exist in our triple- 
parabolic model. Therefore, the extent of these anoma- 
lous divergences of the radius and the interfacial width of 
the critical nucleus depend strongly on the model used. 
Of course, such an anomaly will not be observed in a 
real experiment as the meaning of nucleation becomes 
already obscure near the spinodal as the work of forma- 
tion is of the order of the thermal energy ~ /csT near 
the spinodaP". 

FigureHfa) shows the evolution of the free energy from 
the CNT regime near the two phase coexistence to the 
spinodal regime for the case-Ill. The corresponding evo- 
lution of the density profile is shown in Fig. [DJb). Al- 
though the density profile seems to become flat and dif- 
fuse as the spinodal is approached, in fact the shape does 
not change even though the density difference between 
the inside and the outside of the bubble decreases as the 
positions of the dividing radius as well as the match- 
ing radius rs do not change in this spinodal regime. 

The equimolar dividing radius as the function of the 
scaled undersaturation A^/A/^spin is shown in Fig. El^a) 
and the scaled radius r/j/rmin are shown in Fig. EJb), 
where r^in is the minimum of which is the values at 
the spinodal regime which is in fact the constant. In con- 
trast to the numerical results of Shen and Debenedetti^i 
for the Lennard-Jones fluid, the scaling relation for the 
dividing radius rjj is barely satisfied in the CNT regime 
for the present triple-parabolic free energy model as three 
curves do not fall into a single curve as shown in Fig.lSfb). 
The scaling, however, recovered in the spinodal regime as 
rn is constant. 
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FIG. 4: (a) The evolution of the triple-parabolic free en- 
ergy from CNT regime to the spinodal regime for the case- 
III. (b) The corresponding evolution of the density profile of 
the critical bubble. It seems that the interface of the bubble 
changes from sharp interface to diffuse one as the spinodal is 
approached. 

Shen and Debenedetti^i also noticed that the density 
at the origin 0org 

0org = ^ 0) (29) 

as the function of scaled undersaturation A/i/A/igpin is 
represented by a single curve, and satisfies the scal- 
ing. Furthermore, the curve shows shallow minimum at 
the CNT-spinodal boundary. They^i argued that this 
anomaly is due to the change of the nucleation behavior 
from CNT-nucleation to the spinodal nucleation. Similar 
anomalies are also found for various radii and material 
excess of the critical nucleus in the self-consistent-field 
(SCF) model of binary polymer blendfSS, 

They21 also found that the mean density 0mcan which 
is defined as the average density calculated by integrating 
the bubble density (f)(r) from the origin to the equimolar 
dividing radius td 

0mcatt= / 0(r)Wdr/— ^ (30) 
Jo -J 

also shows the scaling rule, and is expressed by a sin- 
gle curve as the function of the scaled undersaturation 
A/z/A/^spi„. 
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anomalous minimum found by them^ is not due to the 
change of the style of nucleation to anomalous spinodal 
nucleation but merely due to the characteristics of the 
Lennard-Jones fluid. 




FIG. 7: The mean density </> 
undersaturation A/x/A/Xspin. 



as the function of the scaled 



FIG. 5: (a) The equimolar dividing radius td as the functions 
of the scaled undersaturation A/.i/A/.ispin for the three cases, 
(b) The test of the scaling relation for the ro/r-mi-n- The three 
curves do not overlap in the CNT regime. Scaling is recovered 
only in the spinodal regime. 




FIG. 6: The density at the origin (/>org as the functions of 
the scaled undersaturation A/x/A/.ispin for case-I to case-Ill. 
Note that the bulk liquid density \s (j)2 = 1 and the bulk vapor 
density is 0o = 0.1 (Table HI. 



Figure 6 shows the density at the origion (porg of our 
square-gradient model for the three cases. In contrast to 
the results of Shen and Debenedetti^i the scaling rela- 
tion is marginally satisfied. In particular, we cannot find 
out the anomalous density minimum found previously^^. 
In fact, since both calculations of ours and of Shen and 
Debenedetti is based on the density functional theory 
and is limited within the mean-field approximation, the 
anomalous spinodal nucleation is beyond the scope of 
our study as well as of Shen-Debenedetti's. Thus the 



Figure [7] shows the mean density (/)mcan of the bubble 
as the function of the scaled undersaturation A/i/A/igpin- 
Again, the scaling relation is marginally satisfied in the 
CNT regime in contrast to the numerical results by Shen 
and Debenedetti2i. However, the scaling is recovered in 
the spinodal regime, which is due to the facts that the 
equimolar dividing radius rjj becomes constant and the 
density profile becomes universal curve in the spinodal 
regime which will be discussed later in this section. 

It has been argued for many year o^i^^'^^i'^° that the 
interfacial thickness of the critical nucleus diverges as 
the spinodal is approached. Therefore, the nucleation 
becomes anomalous near the spinodal as the critical nu- 
cleus becomes ramified fractal objcct^^ rather than the 
compact spherical droplet. Shen and Debenedetti^i cal- 
culated the interfacial thickness of the critical bubble 
of the Lennard-Jones fluid defined by the 10-90 width 
w^go-io = ''go - '^lo through 

h - (02 - 0org) X 0.1 = (/) (rgo) (31) 

where i^org is defined in Ea. (|29p . and 

'/>org + (02 - 0org) X 0.1 = (rio) . (32) 

They found that the thickness wgo-io of the Lennard- 
Joes fluid does diverge as the spinodal is approached. 
Furthermore, the thickness of the Lennard-Jones fluid 
exhibits scaling. By dividing minimum thickness, the 
interfacial thickness is represented by a single universal 
curve if it is plotted against the scaled undersaturation 

A^/A/igpin. 

We showed the interfacial thickness of our square- 
gradient density-functional model derived from wgo-io — 
^90 ~ ^10 in Fig. [51^ a) as the function of the scaled un- 
dersaturation A/i/A/ispin. By careful examination of 
Eqs. dnj) and with Eqs. ((HI) and ([13, we can easily 
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FIG. 8: (a) The interfacial thickness calculated from 
u)90-io = ''go — ''10 as the functions of the scaled undersat- 
uration A/x/A/^spin. (b) The scaled thickness uigo-io/uispin 
divided by the thickness in the spinodal regime. 



FIG. 9: (a) The density profile </>(r) of the critical bubble for 
the several values of the undersaturation Ap in the spinodal 
regime for Case-II. (b) The corresponding scaled density pro- 
file (j){r) defined by Eq. (|33|l . All curves in (a) overlap into a 
single universal curve. 



prove that both the radius rgo and rio are constant in 
the spinodal regime. Then, so is the interfacial thickness 
Wgo- lo- 
in fact, the thickness wgo_io = ''90 ~ ''10 becomes con- 
stant in the spinodal regime but the magnitude differs 
from case-I to case-Ill in Fig. |Sl[a) . The scaling relation 
is also tested in Fig.[8Ub) by dividing the interfacial thick- 
ness W90-10 by the width Wgpin in the spinodal regime. 
The three distinct curves in Fig.[5]Jb) show that the scal- 
ing for the interfacial thickness is strongly violated in 
the CNT regime in contrast to the results of Shen and 
Debenedetti^i for the Lennard-Jones fluid. The scaling, 
however, is recovered in the spinodal regime as the thick- 
ness W90-10 = Wspin becomes constant. 

In Fig. mja) we show the density profiles of critical 
bubble for several values of undersaturation A/i in the 
spinodal regime. The density profile changes continu- 
ously and it becomes flat and diffuse as the spinodal is 
approached as the density difference between the inside 
and the outside of the bubble becomes small. The fact is 
that the interfacial thickness is constant as the radius rgo 
and rio are constant in this regime, which can be further 
confirmed by inspecting the scaled density profile 0(r) 
defined by 



i^org 



(33) 



<P2 - (Poi-g 

which can be shown to become a universal curve inde- 



pendent of undersaturation A/i in the spinodal regime: 



, , , , /sin (Fir) 
(r) = (/)oA2rBCSc(rirB) I — ^ - Fi 



1 A i\ I exp {-T2r + T2rB) 
l-</>o|Ai|rB (34) 



where 



I Ai I + A2 - A2rBCSC (Firs) Fi 



(35) 



which does not change the shape as the matching radius 
rs does not depend on the undersaturation A/i in the 
spinodal regime. 

We show in Fig. [9{b) the universal scaled density pro- 
file in the spinodal regime. It is clear from this figure 
that the interfacial width is constant up to the spinodal 
in the spinodal regime. 

Cahn and Hilliard^^ used the radius ri/2 defined by 



-f 02)/2 = (/.(ri/2) 



(36) 



as the typical size of the critical nucleus, and has found 
that ri/2 diverges as the spinodal approached. Similar di- 
vergence for 7-1/2 was found for other models^i^. How- 
ever, it is clear from Fig. [H^b) that ri/2 of our triple- 
parabolic feee energy model does not diverge but remains 
constant similar to vb and td in the spinodal regime. 
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FIG. 10: (a) The reversible work of formation W* of the 
critical bubble compared with the classical work of formation 
Wc\ as the functions of the scaled undersaturation A/x/A/ispin. 
(b) The scaled work of formation W,/Wci which does not 
exhibit the scaling property. 



Therefore, the diverging critical nucleus and the di- 
verging interfacial thickness and corresponding fractal 
critical nucleus predicted originally by Cahn and Hilli- 
rad2^ and confirmed by other o^^i^^'^^ cannot be applied 
to the critical bubble of our model. Using the nucleation 
theorem^^i^, Wilcmski and Li^ showed that the excess 
number of molecule 



Ag ~ 47r / [(l>2 — (pir)] r'^dr, 



(37) 



which is in fact the missing number of molecule in the 
bubble nucleation, is proportional to the isothermal com- 
pressibility K of the metastable liquid 



Ag oc K. 



(38) 



Since the isothermal compressibility k diverges as the 
spinodal is approached in the square-gradient density 
functional model with the polynomial free energy con- 
taining up to (j)'^ {(p'^-tield model) ^^'^^ or the Lennard- 
Jones liquid^ii^, the corresponding Ag diverges and so 
does the size of critical nucleus. Now, in our triple- 
parabolic free energy model, the isothermal compressibil- 
ity given by = 0|A2 stays constant as the spinodal 
is approached, which is the reason why the size of bub- 
ble represented by the width wgg-io of the interface does 
not diverge. Therefore our triple-parabolic free-energy 
model that is similar but is slightly different from 0*- 



field model does not give the diverging interfacial thick- 
ness as the spinodal is approached. Therefore, this type 
of parabolic modet^ i^°'^^ may not be suitable to study the 
structure of nucleus near the spinodal. This divergence at 
the spinodal is caused by the divergence of the isothermal 
compressibility of the metastable liquid as the spinodal 
is approachedi^. However, the degree of this divergence 
should depend strongly on the mathematical structure 
of the model used, or the physical system experimentally 
studied. In fact, in contrast to the prediction of Cahn 
and Hilliard^^, the experimental results by Lefebvre et 
al.'^^ for polymer mixtures indicates that the size of crit- 
ical nucleus monotonically decrease rather than diverges 
with increasing quench depth, and it remains finite at the 
spinodal. 

Finally, Fig. [TUf a) shows the reversible work of for- 
mation of the critical bubble calculated from Eqs. ([^0]) 
and ((22l) compared with the classical work of formation 
Ea. (j24p . Our square-gradient density- functional model 
correctly predicts that the work of formation vanishes at 
the spinodal while the classical work remains finite at the 
spinodal. The scaled work of formation W* /Wd barely 
shows scaling properties 



= 1 - 



spin 



(39) 



proposed by Talanquer-i? and Shcn and Debenedeti^i in 
contrast to the results for the Lennard-Jones fluid^^-. The 
shape changes from convex up for the Case-I, which is 
similar to previous worki^i^i for the Lennard-Jones and 
the square-well fluids, to the convex down for the Case- 
III, which is simlar to the results by Miiller et alM- for 
the density fucntional culculation of bubble nucleation 
in the polymer mixture. Therefore, the various scaling 
rules^i^i^i^i^ for the work of formation could depend 
strongly on the model the authors has chosen. Then the 
scaling rule, such as Eq. ((39)) may not always be valid 
for any model or any physical systems. A similar modest 
scaling to ours was found by Kusaka^**. 

Although we do not obtain anomalous structural prop- 
erties of the critical bubble near the spinodal, the work 
of formation of the critical bubble correctly vanishes 
at the spinodal. Therefore anomalous spinodal nucle- 
ationii^i^ should occur even within our model as soon 
as the work of formation is of the order of the thermal 
energy W* ~ ksT. 



IV. CONCLUSION 

In this study, we use a square-gradient density- 
functional model with triple-parabolic free energy, to 
study the various properties of the critical bubble of ho- 
mogeneous bubble nucleation. In particular, we pay at- 
tention to the scaling properties of the reversible work 
of formation of critical bubble as well as other quanti- 
ties such as the density at the center of critical bubble 
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or the interfacial thickness. In contrast to the similar 
density-functional study of Shen and Debenedetti^^ who 
have found that not only the work of formation but other 
various quantities related to the critical bubble shows 
scaling rule, we found marginal scaling relation for our 
model. Since our model is specified only by the com- 
pressibility and the undersaturation, our result seems to 
indicate that the liquid compressibility changes concomi- 
tantly with the vapor compressibility as the temperature 
is changed in the Lennard- Jones fluid, which would lead 
to the scaling rule found by Shen and Debenedetti^^ for 
different temperatures. 

Furthermore, we found that the structural anomalies 
of the critical bubble near the spinodal found in the 
(^''-field theory^^'^° does not exist. Specifically we did 
not find anomalous minimum in the density at the cen- 
ter of critical bubble found by Shen and Debenedetti^i 
for the Lennard-Jones fluid when the undersaturation is 
increased toward the spinodal. They^^ attributed this 
anomaly due to the spinodal nucleation. We also found 
that the interfacial thickness of the critical bubble in the 
spinodal regime does not diverge near the spinodal as 
Unger and Klein^S predicted. In contrast, the thickness 
remains constant up to the spinodal. Furthermore, the 
interfacial profile remains the same in the spinodal regime 



as we scale the profile. Therefore a part of the anomalous 
behavior of various quantities near the spinodal predicted 
by various authors using various different models could be 
largely model-dependent. Such anomalies and their ex- 
tent could partly be due to the mathematical structure of 
(j)*-tie\d or the Lennard-Jones fluid. Further theoretical 
study to test the scaling properties of the critical nucleus 
and bubble, in particular near the spinodal using other 
realistic models will be certainly necessary. The criti- 
cal cavity^"^'-- in the stretched fluid can also be studied 
using this square-gradient density-functional model with 
triple-parabolic free energy. Finally, the dynamics of the 
bubble formation will also be interesting, for that pur- 
pose numerically efficient cell dynamics systent^ could 
be used. 
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